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A  new  theoretical  model  for  the  convergent  nozzle  ejector  in  the  anode  recirculation  line  of  the  polymer 
electrolyte  membrane  (PEM)  fuel  cell  system  is  established  in  this  paper.  A  velocity  function  for  analyzing 
the  flow  characteristics  of  the  PEM  ejector  is  proposed  by  employing  a  two-dimensional  (2D)  concave 
exponential  curve.  This  treatment  of  velocity  is  an  improvement  compared  to  the  conventional  ID  “con¬ 
stant  area  mixing”  or  “constant  pressure  mixing”  ejector  theories.  The  computational  fluid  dynamics 
(CFD)  technique  together  with  the  data  regression  and  parameter  identification  methods  are  applied 
in  the  determination  of  the  velocity  function’s  exponent.  Based  on  the  model,  the  anode  recirculation 
performances  of  a  hybrid  PEM  system  are  studied  under  various  stack  currents.  Results  show  that  the 
model  is  capable  of  evaluating  the  performance  of  ejector  in  both  the  critical  mode  and  subcritical 
mode. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  polymer  electrolyte  membrane  (PEM)  fuel  cell  is  consid¬ 
ered  the  best  candidate  to  replace  the  combustion  power  generator 
because  of  its  capability  of  high  power  density,  low  operating  tem¬ 
perature,  and  short  start-up  time.  The  fuel,  hydrogen,  is  usually 
stored  in  a  highly  pressurized  tank  in  order  to  increase  its  volume 
density,  while  the  pressure  of  the  fuel  cell  stack  is  relatively  low.  The 
high-pressure  difference  between  hydrogen  tank  and  fuel  cell  stack, 
which  contains  abundant  pressure  potential  energy,  can  be  utilized 
by  means  of  an  ejector.  With  no  moving  parts  and  less  maintenance, 
the  ejector  uses  the  high-pressure  hydrogen  as  the  primary  fluid  to 
suck  the  anodic  exhaust  (secondary  flow).  Unconsumed  hydrogen 
in  the  anodic  exhaust  is  recycled  to  the  fuel  cell  stack  so  that  the 
fuel  utilization  efficiency  is  increased  [1-3]. 

Unlike  ejectors  applied  in  refrigeration  and  SOFC  anode  recircu¬ 
lation  systems,  the  ejector  in  the  anode  recirculation  line  of  PEM 
fuel  cell  system  usually  has  a  convergent  not  convergent-divergent 
nozzle  to  convert  the  primary  flow’s  pressure  energy  to  kinetic 
energy.  The  usage  of  the  convergent  nozzle  ejector  is  to  decrease 
the  condensation  of  water  vapor  in  the  ejector  due  to  the  low  tem¬ 
perature  of  the  primary  and  secondary  flows  in  the  PEM  fuel  cell 
system.  A  clear  understanding  of  the  working  principle  and  perfor¬ 
mance  characteristic  of  a  convergent  nozzle  ejector  is  required  for 
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the  performance  evaluation  and  analysis  of  the  anode  recirculation 
PEM  fuel  cell  system. 

The  ejectors  are  widely  used  in  industry  and  have  been  stud¬ 
ied  for  decades  using  theoretical,  empirical  or  numerical  methods 
[4-6].  However,  most  of  the  existing  ejector  models  are  developed 
based  on  one-dimensional  (ID)  technique  for  cooling  and  refrig¬ 
eration  [7-9]  and  SOFC  [10-13]  applications.  There  are  very  few 
researches  for  the  ejector  applied  in  the  fuel  recirculation  line  of 
the  PEM  fuel  cell  system. 

Karnik  et  al.  studied  an  ejector  based  PEM  fuel  cell  anode 
recirculation  system.  The  ejector  was  modeled  based  on  the  ID  sim¬ 
plification  and  the  “constant  area  mixing”  assumption  [14,15].  Bao 
et  al.  proposed  models  for  a  PEM  fuel  cell  system  with  an  ejector 
based  recirculation  line,  however  only  the  ejector  performances  in 
the  critical  operational  condition  were  analyzed  [  16].  He  et  al.  inves¬ 
tigated  a  hybrid  fuel  delivery  system  in  the  PEM  fuel  cell  consisting 
of  two  supply  and  two  recirculation  lines.  The  ejector  was  assumed 
work  in  the  critical  mode  and  the  primary  and  secondary  flows 
were  analyzed  based  on  the  ID  “constant  pressure  mixing”  theory 
[17].  According  to  the  ejector  flow  theory,  the  ejector  performance 
is  divided  into  three  operational  modes,  i.e.,  back  flow,  subcritical 
and  critical  modes  by  the  operating  conditions  [4,8].  The  ejector 
may  work  in  the  subcritical  mode  or  even  back  flow  mode  during 
start  up,  load  changes  and  shut  down.  In  these  cases,  the  flow  char¬ 
acteristic  becomes  complicated  and  unexpected  fluctuations  in  the 
anode  recirculation  line  could  occur. 

In  this  paper,  we  aim  to  develop  a  simple  and  accurate  theoret¬ 
ical  model  for  both  critical  and  subcritical  operating  ejector  in  the 
hybrid  PEM  fuel  cell  system.  The  paper  is  organized  as  follows: 
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Nomenclature 

A  area(m1 2) 

Cp  specific  heat  of  gas  at  constant  pressure  (J  kg-1  K-1 ) 
D  diameter  (m) 

m  mass  (kg) 

M  Mach  number 

Mo  molecular  weight  (kg  mol-1 ) 

n  molar  flow  rate  (mol  s-1 ) 

nv  exponent  of  the  2D  curve 
P  pressure  (Pa) 

r,  R  radius  (m) 

Rg  gas  constant  (J  kg-1  K-1 ) 

Ru  universal  gas  constant  (J  mol-1  K-1 ) 

T  temperature  (K) 

v,  V  velocity  (ms-1 ) 

W  mass  flow  rate  (kg  s-1 ) 

y  mass  fraction 

Greek  symbols 

1 6P  pressure  ratio,  (pressure  in  bar) 

Pd  diameter  ratio,  Dm/Dt 

k  specific  heat  ratio  of  gas 

co  recirculation  ratio,  WS/WP 

p  density  (kg  nrr3) 

*Pp  isentropic  coefficient  of  primary  flow 

&S  isentropic  coefficient  of  secondary  flow 
§exp  coefficient  accounting  for  friction  loss  during  the 
mixing  process 

Subscripts 

an  anode  flow  channel 

MEA  membrane 

P  primary  flow  (i.e.  inlet  fuel) 

S  secondary  flow  (i.e.  anode  recycling  gas) 

B  ejector  exit 

t  nozzle  throat 

0  ejector  inlet 

1  nozzle  exit 

2  mixing  chamber  inlet 

3  ejector  exit 

in  inlet 

out  outlet 

Superscript 

i  chemical  component 


1.  The  working  principle  of  an  ejector  based  anode  recirculation 
PEM  fuel  cell  system  is  described  and  the  PEM  fuel  cell  model  is 
briefly  introduced  in  Section  2. 

2.  In  Section  3,  a  new  theoretical  model  for  the  convergent  nozzle 
ejector  is  constructed  by  employing  a  2D  exponential  function 
to  compute  fluid  velocity  near  the  ejector  inner  walls.  Govern¬ 
ing  equations  for  computing  the  mass  flow  rate,  recirculation 
ratio  and  the  ejector  exit  temperature  are  derived  based  on  the 
thermodynamic  and  fluid  dynamic  principles. 

3.  A  method  for  determination  of  the  2D  exponential  function  is 
proposed  in  Section  4.  A  set  of  ejector  geometries  and  work¬ 
ing  conditions  are  first  simulated  using  the  computational  fluid 
dynamics  (CFD)  technique.  Then  the  simulation  results  are 
treated  using  the  data  regression  and  parameter  identification 
methods  to  find  out  the  relation  between  the  exponent  and  the 
geometries  and  working  conditions. 


4.  In  Section  5,  the  detailed  application  procedure  for  the  ejector 
model  is  presented.  Comparisons  are  made  between  the  model 
results  and  the  corresponding  CFD  data. 

5.  In  Section  6,  a  hybrid  PEM  fuel  cell  system  integrated  with 
an  ejector  at  the  anode  recirculation  line  is  studied  using  the 
proposed  model.  Behaviors  of  the  anode  recirculation  line  are 
obtained  and  analyzed. 

2.  PEM  fuel  cell  system 

2.1  Working  principle  of  anode  recirculation  PEM  fuel  cell  system 

A  PEM  fuel  cell  system  with  an  ejector  in  the  anode  recircula¬ 
tion  line  is  shown  in  Fig.  1.  It  mainly  consists  of  three  components: 
a  PEM  fuel  cell  stack,  an  ejector  and  a  humidifier.  The  structure  of  a 
PEM  fuel  cell  can  be  described  as  two  electrodes  (anode  and  cath¬ 
ode)  separated  by  a  membrane.  High-pressure  hydrogen  entering 
the  ejector  is  called  the  primary  flow.  Inside  the  ejector,  the  primary 
flow  expands  in  a  convergent  nozzle  in  which  its  pressure  poten¬ 
tial  energy  converts  into  kinetic  energy,  resulting  in  a  low-pressure 
region.  The  anode  recycle  gas  (secondary  flow)  is  entrained  into 
the  ejector.  The  two  flows  mix  in  the  ejector  then  discharge  to  the 
humidifier.  According  to  Refs.  [14-17],  the  functions  of  ejector  in 
the  anode  recirculation  line  are 

•  Utilize  the  pressure  potential  energy  of  hydrogen  otherwise 
wasted; 

•  Recycle  the  unconsumed  hydrogen  to  the  fuel  cell  to  increase  the 
fuel  usage  efficiency; 

•  Regulate  the  anode  humidity  with  the  recycle  gas; 

•  Raise  the  secondary  flow  pressure  to  meet  the  fuel  cell  pressure 
at  the  required  level. 

The  ejector  performance  is  usually  evaluated  by  the  recirculation 
ratio  co ,  which  is  defined  as 


co  =  Ws/WP 


(1) 


where  Wp  and  Ws  are  the  mass  flow  rates  of  the  primary  and  sec¬ 
ondary  flows,  respectively. 

In  the  humidifier,  the  gas  humidification  is  considered  to  be  per¬ 
fect.  For  a  volume  of  hydrogen  gas  with  water  vapor  and  liquid 
water,  the  water  activity  aw  and  relative  humidity  are  defined 
as 


mH2o  .  ,  _  Pv 

HI y, sat’  _  Pv,sat 


(2) 


where  mH2o>  rnv.sat  is  the  total  water  and  saturation  mass  of  water 
vapor  in  the  volume,  respectively.  Pv  is  the  partial  pressure  of  water 
vapor  in  the  mixture.  Pv>Sat  is  the  saturation  pressure  of  the  water 
vapor,  which  is  related  to  the  mixture  temperature  [18]. 

In  the  anode,  Hydrogen  dissociates  into  protons  that  flow 
through  the  membrane  to  the  cathode.  At  the  same  time,  electrons 
lost  from  hydrogen  which  are  collected  as  electrical  current  by  an 
external  circuit  linking  the  two  electrodes.  Air  flows  to  the  cathode 
where  oxygen  combines  with  the  electrons.  The  chemical  reactions 
taking  place  at  the  two  electrodes  are  listed  as  follows: 


Anodereaction  :  2H2  ->►  4H+  +4e 

Cathodereaction  :  02+4H++4e  ->  2H20 

Overallreaction  :  02  +  2H2^  2H20  +  heat 


(3) 

(4) 

(5) 


2.2.  PEM  fuel  cell  model 

The  model  for  the  PEM  fuel  cell  consists  of  two  parts:  one  is 
about  thermodynamic  equations  to  calculate  the  temperature  and 
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mass  flow  rate,  and  the  other  is  electrochemical  part  that  computes 
the  cell  potential  and  power.  As  the  main  purpose  of  this  work  is 
to  study  the  performance  of  fuel  delivery  and  anode  recirculation 
lines,  only  the  performance  of  anode  side  is  analyzed.  The  steady- 
state  PEM  fuel  cell  model  adopted  in  this  work  is  briefly  introduced 
as  follows: 

The  mass  conservation  equations  of  the  anode  are 


^H2,an,in  —  ^H2,an,out  +  ^H2,  reacted 

(6) 

Wv,an,in  —  Wv,an,out  —  Wv,MEA 

(7) 

where  WHl  an  in  and  Wv  an  in  are  the  mass  flow  rates  of  hydrogen  and 
water  vapor  entering  the  anode  flow  channel,  which  are  related  to 
the  ejector  and  humidifier’s  operational  conditions: 

^H2,an,in  =  +yH2,an,out^S 

(8) 

i/i/  yv,hum  .. j 

''M/,an,in  —  *M-I2,an,in 

1  -yv,hum 

0) 

where  WP  and  Ws  in  Eq.  (8)  are  derived  from  the  ejector  model 
presented  in  Sections  3-5.  yv,hum  is  the  mass  fraction  of  the  water 
vapor  in  the  humidifier. 

In  Eq.  (6),  WH2  reacted,  which  means  the  mass  rate  of  hydrogen 
consumed  in  the  anode,  is  described  as  [19]: 

/Moh-, 

Wh2, reacted  =^-2^  (10) 

WvMEA  in  Eq.  (7)  is  the  transport  rate  of  water  vapor  through  the 
membrane  from  the  anode  to  the  cathode  resulting  from  diffusion, 
electro-osmotic  drag  force,  and  pressure  gradient: 

/Moh-,0 

WvMEA  =  anetNcen—p^-  (11) 

where  anet  is  the  net  water  transfer  coefficient  which  is  defined  as 
the  ratio  of  the  net  water  flux  in  PEM  to  the  ion  flux  [19,20]. 

The  expression  of  the  voltage  of  a  single  cell  is 

Vcell  =  E  +  T)ac  t  +  TJ  ohmic  +  Pdiff  (12) 


Pdiff  =  m  exp(n/) 


(16) 


where  the  parameters  §l-§4,  m  and  n  can  be  found  in  literatures 
[21-24]. 

The  thermodynamic  efficiency  of  the  PEM  fuel  cell  can  be  eval¬ 
uated  by  [25]: 


Eff  = 


2nce\\Vce\\I 

Wh2,  reacted  '  LEIVh2 


(17) 


3.  Ejector  theoretical  model  development 

A  convergent  nozzle  ejector  applied  in  the  anode  recirculation 
line  of  PEM  fuel  cell  is  schematically  shown  in  Fig.  2(a).  The  ejec¬ 
tor  has  a  convergent  nozzle  not  convergent-divergent  nozzle  as 
the  flow  expansion  device.  According  to  the  characteristic  of  the 
secondary  flow,  the  ejector  performance  are  divided  into  three 
operational  modes:  back  flow,  subcritical  and  critical  as  depicted 
in  Fig.  2(b)  [4,8].  The  secondary  mass  flow  rate,  which  is  very  sensi¬ 
tive  to  the  operational  conditions  in  the  subcritical  mode  and  keeps 
near  constant  in  the  critical  mode. 

Without  loss  of  generality,  the  following  assumptions  are  made 
in  developing  a  model  for  the  convergent  nozzle  ejector  used  in  the 
PEM  fuel  cell  system: 


1.  The  primary  flow  is  treated  as  the  ideal  gas. 

2.  The  primary  flow  velocity  is  uniform  in  the  radial  direction  inside 
the  ejector. 

3.  The  velocity  of  the  secondary  flow  inside  the  ejector  is  non- 
uniformly  distributed  in  the  radial  direction  and  there  exists 
velocity  boundary  layer  near  ejector  inner  walls. 

4.  Pressure  and  temperature  of  both  the  primary  and  the  secondary 
flows  are  uniformly  distributed  in  the  radial  direction  of  ejector. 

5.  The  isentropic  relations  hold  in  calculating  friction  losses. 


3.1.  Primary  flow 


where  the  thermodynamic  potential  of  the  H2  +  02  reaction,  E,  can 
be  further  expressed  by  [21  ]: 

E  =  1.229  -  0.85  x  10“3(T  -  298.15) 

+  4.3085  x  10-5T  ln[p*2(p*2)0-5]  (13) 

The  parametric  expression  of  the  activation  over  potential  pact, 
ohmic  over  potential  pohmic  and  diffusion  over  potential  pdiff  can  be 
written  as  follows: 

(14) 

(15) 


According  the  flow  characteristic  of  the  convergent  nozzle,  the 
flow  through  the  convergent  nozzle  is  divided  into  two  different 
regions:  subsonic  and  sonic  flow  by  the  pressure  ratio  Ps/Pp.  The 
relation  between  the  pressure  ratio  and  the  mass  flow  rate  is  shown 
in  Fig.  3.  Due  to  without  a  divergent  part,  the  flow  cannot  reach  sonic 
flow  condition  when  the  pressure  ratio  is  greater  than  the  critical 
value  ucr.  The  critical  pressure  ratio  ucr  of  a  convergent  nozzle  can 
be  written  as 

Utt)  (1S) 

For  sonic  flow  (when  PP0  >Ps,o/DCr),  the  primary  flow  rate  mP 
can  be  obtained  by  using  the  isentropic  flow  relations  and  energy 


Pact  =  §1  +  §2^11  +  ^’f’cell  ln(Q)2)  +  ?4^cell  in(0 

~  j  pinternal 

Pohmic  =1-  K 
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(a)  Suction  chamber 

Mixing  chamber 

Diffuser 

Secondary  flow 

_ 

Primary  flow  ■ - 

:A 

r.  A  D 

v  m 

- /  - r 

=  Convergent  nozzle  : 

0  12  3 


Fig.  2.  Schematic  diagram  showing  geometries  and  operational  modes  of  an  ejector  (a)  Ejector  geometries;  (b)  Operational  modes. 


balance  law: 


v  (k+1)/(2(k-1)) 


Wp=PpAt\l^rP(-K+-u 

and  the  Mach  number  at  the  nozzle  throat  is  1,  i.e.: 


For  subsonic  flow  (when  Pp  <PslvCr),  the  primary  flow  rate  mp 
and  the  Mach  number  at  the  nozzle  throat  Mt  are  calculated  from: 


(19) 


Wp  =  PpA±\ 


1 2ilrPK[(Ps/Pp)2lK  -  (Ps/Pp)(1+'0A] 


(K-l)RgTp 


(21) 


Mt  =  1 


(20) 


(22) 

where  *Pp  is  the  isentropic  coefficient  taking  into  account  the  flow 
friction  loss. 

The  primary  flow  expands  fully  in  the  suction  chamber.  The 
ambient  pressure  of  the  expansion  flow  can  be  represented  by  the 
secondary  flow  pressure  Ps.  Using  the  isentropic  flow  and  energy 
balance  laws  for  the  primary  flow  from  Section  1  to  Section  2,  we 
have: 


Pp 

Ps 


1  +  2  (*■  —  1  )Mp2 


W(K-l) 


2^  =  l  +  l(/f-l)M,2 

Ip, 2  2. 

Vp2  =  MP2y^cRgf^2 


P,2 


§expUp,2 

Dt 


2  +  (/c-l)Mp22l(-1)/(4(^1)Y  Mt 

2  +  (/c-l)M2J  \MPa 


0.5 


(23) 

(24) 

(25) 

(26) 


where  §exp  is  a  coefficient  accounting  for  the  frictional  loss  due  to 
the  mixing  process  of  the  two  flows. 
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3.2.  Secondary  flow 


in  Section  2 

2Vp?2(i?2  +  ftp, 2  +  rcvftp,2) 


Vs,2  = 


^s,2(nv  +  l)(nv  +  2)(i?2  +  Kp, 2) 
where  the  average  density  of  the  secondary  flow  is  given  as 


Ps  = 


Ps  Ps 


J2nsMoi 


(31) 


(32) 


3.3.  Energy  balance  of  two  flows 

For  an  ideal  gas,  the  energy  balance  of  the  primary  and  secondary 
flow  in  the  ejector  can  be  described  by 

EWp  +  JfwlCpTs  =  ]T(W'  +  W's)CpT3  +  Eloss  (33) 

i  i  i 


The  ejector  performance  is  significantly  affected  by  the  flow 
characteristics  in  Section  2  [9].  In  Section  2,  a  mixing  layer  sepa¬ 
rates  the  primary  flow  and  the  secondary  flow.  The  secondary  flow 
outside  the  layer  has  a  non-linear  velocity  distribution  due  to  the 
turbulent  flow  and  fluid  viscosity  as  shown  in  Fig.  4.  In  the  conven¬ 
tional  1 D  analysis,  both  velocities  of  the  primary  flow  and  secondary 
flow  are  treated  uniform  in  the  radial  direction.  Obviously,  this 
method  will  induce  uncertainty  in  the  ejector  modeling. 

It  is  important  to  note  that  the  velocity  distribution  in 
the  convergent  nozzle  ejector  is  different  from  that  in  a 
convergent-divergent  nozzle  ejector.  The  velocity  of  the  primary 
flow  at  nozzle  exit  is  relatively  low  in  the  convergent  nozzle  ejec¬ 
tor  due  to  no  divergent  part.  Therefore,  the  primary  flow  in  the 
convergent  nozzle  ejector  cannot  accelerate  the  secondary  flow 
effectively  like  the  convergent-divergent  nozzle  ejector  do,  leading 
to  a  concave  velocity  distribution  of  the  secondary  flow  as  depicted 
in  Fig.  4. 

In  order  to  describe  the  concave  velocity  distribution  curve,  we 
propose  the  following  governing  equations  for  the  velocities  of  pri¬ 
mary  flow  and  secondary  flow  in  Section  2 


where  the  energy  loss,  Eioss,  of  the  primary  and  secondary  flows  in 
the  ejector  can  be  approximated  as 

floss  =  \o  -  ft>)WpV£2  +  1(1  -  fsWsVf  (34) 

Remark.  The  velocity  of  the  secondary  flow  is  accurately  mod¬ 
eled  by  a  2D  exponential  function  in  the  new  model.  Because  of 
this  improvement,  the  present  model  is  capable  of  predicting  the 
ejector  performance  within  less  uncertainties  compared  to  the  con¬ 
ventional  ID  analysis.  Additionally,  the  treatment  of  the  secondary 
flow  velocity  makes  the  new  model  simpler  than  the  conventional 
“constant  area  mixing”  or  “constant  pressure  mixing”  theory.  For 
example,  the  model  for  calculating  Ws  only  consists  of  one  algebraic 
equation  (Eq.  (30)). 

Note  that  the  exponent  nv  in  Eq.  (30)  is  very  important,  which 
decides  the  pattern  of  the  concave  velocity  distribution  curve.  The 
value  of  nv  is  not  a  constant,  which  will  be  influenced  by  the  ejec¬ 
tor  geometries  and  working  conditions.  In  order  to  determine  the 
exponent  nv,  a  numerical  method  based  on  the  CFD  technique  is 
proposed  in  Section  4. 


Vp  2  (0  <  r  <  Rp  2) 

Vp,2  2_  ftp  2  ^  (Pp,2<r<R2) 


(27) 


where  nv  is  the  exponent  of  the  velocity  function.  nv  is  supposed  to 
be  greater  than  1  (nv  >  1 )  for  a  convergent  nozzle  ejector.  The  graph 
of  vr  versus  r  is  given  as  the  dashed  curve  in  Fig.  4.  Clearly,  it  is  closer 
to  the  actual  velocity  distribution  than  the  ID  approach. 

Based  on  this  2D  velocity  function,  we  define  the  mean  mass 
flow  rate  of  secondary  flow  at  Section  2  as 
rR  2 

Ws=  pvrdA  (28) 

7j?p  2 

where  p  stands  for  the  average  density  of  secondary  flow.  Substi¬ 
tuting  Eq.  (27)  into  Eq.  (28)  leads  to  the  following  integral  equation: 

Ws  =  rR27rVn2^nv  f  2  (r2  -  r)"vrdr  (29) 

(r2~rp,2)  JRp2 

By  evaluating  the  integral  of  Eq.  (29),  the  mass  flow  rate  of  the 
secondary  flow  can  finally  be  expressed  as 


Ws 


27TpsVpp{R2  ~  Pp,2)(P2  +  Pp,2  +  nvPp,2) 

(nv  +  1  )(nv  +  2) 


(30) 


By  using  the  conditions  of  Vs,2  =  Ws/(ps71s,2)  and  As>2  =  tt  {R%  - 
Rp2),  we  obtain  the  average  velocity  of  the  secondary  flow 


4.  Velocity  exponent  determination 

In  this  section,  the  numerical  CFD  method  is  first  introduced 
including  the  grid  structure  building,  solution  strategy  construct¬ 
ing,  turbulence  model  selecting  and  boundary  conditions  setting. 
The  CFD  method  for  the  ejector  flow  simulation  has  been  exper¬ 
imentally  validated  previously  [26].  Based  on  this  method,  four 
different  ejector  geometries  and  32  working  conditions  are  sim¬ 
ulated.  Finally,  the  data  regression  and  parameter  identification 
methods  are  used  to  determine  the  velocity  exponent  nv  based  on 
the  obtained  simulation  results. 

4.1.  Numerical  CFD  method 

The  CFD  method  is  capable  of  producing  details  of  the  flow 
field  and  fluid  properties  based  on  numerical  solutions  of  the  flow 
domain.  The  flow,  heat  transfer,  radiation,  turbulence,  etc.  at  any 
given  operating  conditions  and  model  geometries  can  be  simu¬ 
lated.  Data  that  is  difficult  to  obtain  in  an  experiment  can  be  easily 
analyzed  using  CFD. 

The  commercial  software  Gambit  2.2  and  FLUENT  6.2  are  used 
as  the  grid  generator  and  the  CFD  solver  in  this  numerical  study, 
respectively.  The  ejector  geometries  are  modeled  in  a  2D  domain 
in  the  Gambit.  The  grids  are  only  adapted  at  the  locations  with  sig¬ 
nificant  flow  changes  such  as  velocity  boundary  and  shock  position 
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Table  1 

Four  ejector  geometries. 


for  faster  computation  speed.  3D  effects  can  be  reflected  by  the  2D 
ejector  model  since  the  axi-symmetric  solver  is  applied. 

The  flow  inside  the  ejector  is  governed  by  the  compressible 
steady-state  turbulent  form  of  the  flow  governing  equations.  The 
non-linear  governing  equations  are  solved  using  the  “Segregated- 
implicit”  solver  [27].  The  species  transport  model  is  used  to  deal 
with  the  mixture  fluid  flow.  The  second-order  upwind  scheme  is 
adopted  to  discretize  convective  terms.  The  SIMPLEC  algorithm  is 
applied  to  obtain  the  pressure  field.  Moreover,  the  “RNG  k-e  model” 
is  selected  as  the  turbulence  model  for  its  ability  to  better  predict 
ejector  performance  than  other  turbulence  models  [26]. 

The  working  fluid  of  the  primary  flow  is  pure  hydrogen,  while 
the  secondary  flow  is  composed  of  hydrogen  and  water  vapor.  The 
density  is  obtained  using  the  ideal  gas  relationship.  Boundary  con¬ 
ditions  of  the  primary  flow  and  the  secondary  flow  inlets  are  set 
as  “pressure  inlet”  condition,  and  the  “pressure  outlet”  condition  is 
adopted  on  the  ejector  outlet.  For  each  simulation,  the  solution  is 
iterated  until  convergence  is  achieved  (residue  for  each  equation  is 
less  than  10-4). 

4.2.  Data  regression 

The  pressures  PP  and  Ps  and  diameters  Dt  and  Dm  affect  the 
velocity  distribution  inside  the  ejector,  so  as  to  decide  the  velocity 
exponent  nv.  In  order  to  find  out  the  relations  between  the  pres¬ 
sures,  diameters  and  nv,  we  first  investigate  the  influence  of  PP  and 
Ps  on  nv  at  the  same  ejector  geometries,  then  study  the  effect  of  Dt 
and  Dm  at  constant  working  conditions. 

Four  ejectors  geometries  are  created  as  shown  in  Table  1.  In 
these  ejectors,  the  diameter  ratio  /3D(Dm/Dt )  is  in  the  range  of 
2.25-2.86.  According  to  the  typical  operating  conditions  of  PEM  fuel 
cells,  three  pressure  groups  are  designed:  PB  =  3  bar  and  Ps  =  2.8  bar; 
PB  =  1.5  bar  and  Ps  =  1.35  bar;  PB  =  1.1  bar  and  Ps  =  1.0  bar.  For  each 
pressure  group,  the  ejectors  are  simulated  by  varying  the  primary 
flow  pressure  PP.  The  outputs  of  the  CFD  model  simulation  include 
the  two  mass  flow  rates  Ws  and  WP  and  the  ejector  exit  tempera¬ 
ture. 

Substituting  Ws  and  WP  into  Eqs.  (19)-(26)  and  Eqs.  (31)-(32), 
the  corresponding  nv  can  be  computed  by  an  iterative  calculation. 
As  shown  in  Fig.  5,  it  is  found  that  nv  has  a  exponential  relation  with 
Ps°-8/P^  using  the  data  regression  analysis.  Results  in  Fig.  6  show 
that  nv  is  finely  linear  with  Dm/Dt  in  all  the  studied  cases. 

4.3.  Parameter  identification 

Based  on  the  data  regression  studies,  it  is  reasonable  to  assume 
that  nv  has  the  following  expression: 

nv=A!  exp(pP/0.05)  +  A2pD+A3  (35) 

where  pP  =  P°  8/P^A  (here  Ps  and  PP  is  in  bar);  /3D  =  Dm/Dt;  Alt  A2 
and  A3  are  parameters  which  can  be  determined  by  the  identifica¬ 
tion  method  introduced  as  follows: 

If  N  experimental  or  numerical  CFD  tests  are  conducted  for  dif¬ 
ferent  fip,  Pd  and  nv,  Eq.  (35)  falls  into  a  system  of  linear  equations: 


Fig.  5.  Relation  between  nv  and  ejector  working  pressures. 


where  X  is  the  parameters  matrix  to  be  determined,  A  is  the  mea¬ 
surement  or  simulation  error.  Matrixes  ofX,  *P,  r  are  written  as 


x  = 

■*r 

^2 

;  r  = 

;  <p  = 

"  exp(/6j,/0.05) 
exp(/?p/0.05) 

eh  n 
f%  1 

_^3  _ 

_nv  _ 

_exp(^/0.05) 

— i 

(37) 

The  best  estimate  X*  of  X  can  be  found  using  the  standard  least 
squares  method  proposed  by  Young  [28]  as 

x*  =  (^T^)“VTr  (38) 

In  this  work,  based  on  the  four  ejector  geometries  and  32  work¬ 
ing  conditions,  the  parameters  A\ ,  A2  and  A3  are  obtained  using  the 
above  parameter  identification  method.  Substituting  the  parame¬ 
ters  into  Eq.  (35),  the  velocity  exponent  nv  is  finally  achieved  as 

nv  =  1 .393  x  10“4  exp(/3p/0.05)  +  0.456y6D  +  0.1668  (39) 

The  correlation  of  nv  in  Eq.  (39)  is  obtained  based  on  var¬ 
ious  PEM  ejector  operating  conditions  and  geometries  in  the 
ranges  of  1.75  <PP  <  10.5  bar,  1.0  <PS  <2.8  bar,  1.1  <PB  <3.0  bar, 
2.1  <  Dt  <  3.2  mm  and  5.2  <  Dm  <  8.0  mm.  Furthermore,  the  correla¬ 
tion  of  nv  only  consists  of  four  measurable  variables:  Ps,  PP,  Dm  and 


m  =  T  +  A 


(36) 


Fig.  6.  Relation  between  nv  and  ejector  geometries. 
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END 


Fig.  7.  Solution  procedure  of  the  ejector  model. 

Dt.  Therefore,  the  Eq.  (30)  for  calculating  Ws  will  be  convenience  to 
use  in  practice. 

5.  Model  application 

For  a  given  ejector  geometry,  its  performance  depends  on  the 
pressure  PP  and  Ps,  temperature  TP  and  Ts,  chemical  composition  of 
inlet  fuel  and  anodic  recycle  gas.  The  outputs  of  the  ejector  model 
used  in  PEM  fuel  cell  systems  will  be  the  primary  mass  flow  rate  WP, 
the  secondary  mass  flow  rate  Ws,  recirculation  ratio,  temperature 
and  chemical  compositions  at  ejector  exit.  The  solution  procedure 
of  the  ejector  model  is  given  in  Fig.  7. 

Starting  with  the  same  input  conditions  PP,  Ps,  TP,  Ts  andys,H2o, 
results  such  as  WP,  Ws  and  TB  are  obtained  from  the  ejector  model 
and  CFD  simulation.  Comparisons  of  WP,  Ws  and  TB  between  the 
model  and  CFD  simulation  are  presented  Fig.  8(a)-(c),  respec¬ 
tively.  It  is  found  that  the  predicted  WP  and  TB  by  present  model 
match  fairly  well  those  from  the  CFD  simulation  data.  The  average 
derivations  of  WP  and  TB  between  model  and  CFD  simulation  are 
0.63%  and  0.28%,  respectively.  Note  that  the  flow  characteristic  of 
the  secondary  flow  is  complicated  so  that  Ws  is  difficult  to  pre¬ 
dict  particularly  when  the  ejector  works  at  the  subcritical  mode. 
In  Fig.  8(c),  it  can  be  seen  that  the  new  ejector  model  provides 
reliable  results  in  predicting  Ws.  The  average  derivation  between 
model  and  CFD  simulation  is  7.25%  for  all  the  tested  working 
conditions. 

6.  Results  and  discussions 


Fig.  8.  Comparison  of  results  between  CFD  simulation  and  model,  (a)  Mass  flow 
rate  of  the  primary  flow;  (b)  mass  flow  rate  of  the  secondary  flow;  (c)  ejector  exit 
temperature. 


As  an  example,  Table  2  shows  typical  values  for  a  PEM  fuel  cell  ratio  is  chosen  to  be  constant  2,  and  the  relative  humidity  at  the 

system.  The  ejector  1  in  Table  1  is  adopted  in  the  following  analysis.  humidifier  exit  is  kept  at  80%. 

Starting  with  these  data,  the  steady-state  performances  of  the  ejec¬ 
tor  and  fuel  cell  system  are  studied.  In  order  to  simplify  the  analysis,  6.1.  Ejector  performance 
it  is  assumed  that  the  stack  temperature  is  constant  353  K.  All  the 

gases  are  considered  as  the  ideal  gas,  and  the  temperature  variation  Performances  of  the  ejector  and  fuel  cell  stack  can  be  analyzed 

of  gas  properties  is  neglected.  For  simplicity,  the  air  stoichiometric  by  varying  the  stack  current  and  anode  pressure.  In  this  study, 
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Table  2 

Model  parameters. 


Description 

Value 

Cell  active  area  (m2) 

0.1 

Cell  number 

381 

Membrane  thickness  (m) 

127  x 10-6 

Membrane  equivalent  weight  (kg  mol-1 ) 

1.1 

Membrane  density  (kg  m-3 ) 

2000 

Hydrogen  inlet  temperature,  TP  (K) 

298 

Cathode  inlet  oxygen  relative  humidity  (%) 

100 

Anode  inlet  relative  humidity  (%) 

80 

Stack  operating  temperature  (I<) 

353 

Humidifier  temperature  (I<) 

353 

Anode  pressure  (bar) 

3.0 

Cathode  inlet  pressure  (bar) 

3.0 

Anode  inlet  flow  coefficient  (kg  s-1  Pa-1 ) 

0.3  x  10-6 

Anode  outlet  flow  coefficient  (kgs-1  Pa-1 ) 

0.3  x  10-6 

Cathode  inlet  flow  coefficient  (kg  s-1  Pa-1 ) 

2.3  x  10-6 

Cathode  outlet  flow  coefficient  (kgs-1  Pa-1 ) 

2.3  x  10-6 

Volume  of  lumped  cathode  (m3) 

10.8  x  10-3 

Volume  of  lumped  anode  (m3) 

9.6  x  10-3 

eleven  different  stack  currents  are  selected.  Results  for  analysis  of 
the  anode  recirculation  line  are  shown  in  Figs.  9-11. 

The  behavior  of  PP,  Ps  and  PB  at  different  stack  currents  are 
shown  in  Fig.  9.  It  can  be  seen  that  PP  rises  as  the  current  increases, 
evidently  because  raising  PP  will  increases  the  mass  flow  rate  of 


Fig.  9.  Pressures  at  anode  recirculation  line  with  different  stack  currents. 


Fig.  10.  Velocity  exponent  nv  at  different  primary  flow  pressures. 


Fig.  11.  Ejector  performances  at  different  primary  flow  pressures,  (a)  Mass  flow 
rates;  (b)  Ejector  recirculation  ratio  and  exit  temperature. 


the  fuel  hydrogen.  The  behavior  of  PP  for  a  convergent-divergent 
nozzle  ejector  at  the  same  conditions  is  also  depicted  as  the  dash 
line  in  Fig.  9.  Compared  of  PP  between  the  convergent  ejector  and 
convergent-divergent  nozzle  ejector,  it  is  found  that  there  exists  a 
divergent  position  at  about  Ps/Pp  =  0.528.  This  is  in  line  with  the 
thermodynamic  principles  that  the  critical  pressure  ratio  ucr  is 
0.528  for  a  convergent  nozzle  with  a  diatomic  molecule  as  the  work¬ 
ing  fluid.  In  Fig.  9,  the  pressure  lift  in  the  ejector  increases  with  the 
stack  current,  because  the  pressure  loss  in  the  PEM  fuel  cell  has 
direct  ratio  with  the  mass  flow  rate  (i.e.  the  stack  current). 

The  velocity  exponent  nv  is  shown  in  Fig.  10.  It  is  found  that  that 
nv  is  very  sensitive  to  the  primary  flow  pressure  in  the  low  pressure 
regions,  and  becomes  near  constant  at  high  pressure  regions.  This 
is  because  that  the  ejector  works  in  the  subcritical  mode  when  the 
primary  flow  pressure  is  low  as  shown  in  Fig.  2(b). 

The  ejector  performances  at  different  stack  currents  are  shown 
Fig.  11.  The  characteristics  of  WP  and  Ws  in  Fig.  11(a)  are  very 
similar  to  the  previous  experimental  and  numerical  data  in  Refs. 
[6,8].  Wp  in  the  convergent  ejector  is  different  from  that  of  the 
convergent-divergent  nozzle  ejector  when  Ps/Pp  is  greater  than 
0.528.  Ws  remains  quite  constant  in  the  high  primary  flow  pressure 
region. 

The  ejector  recirculation  ratio  oo  and  ejector  exit  temperature  are 
shown  in  Fig.  11(b).  The  ejector  has  a  good  recirculation  capability 
{(l)> 2)  when  the  primary  flow  pressure  is  greater  than  5 bar.  It  is 
important  to  note  that  the  recirculation  performance  will  decline 
once  the  ejector  works  at  subcritical  mode. 
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stack  efficiency  are  illustrated  at  different  currents.  The  stack  effi¬ 
ciency  increases  with  an  increase  of  the  current,  and  the  cell  voltage 
and  power  decrease  as  the  current  increases.  Although  the  ejector 
performances  are  affected  by  the  stack  current  as  in  Fig.  11,  the 
behaviors  of  these  performance  properties  agree  with  the  perfor¬ 
mance  characteristics  of  a  typical  PEM  fuel  cell  stack.  This  is  because 
that  the  hydrogen  stoichiometric  ratio  changes  correspondingly  as 
shown  in  Fig.  13.  From  Figs.  11  and  13,  it  can  be  seen  that  the 
hydrogen  stoichiometric  ratio  has  a  similar  trend  with  the  ejector 
recirculation  ratio.  The  greater  recirculation  ratio,  which  means  the 
more  unconsumed  hydrogen  recycled  into  the  anode,  leads  to  the 
larger  hydrogen  stoichiometric  ratio.  In  order  to  adjust  the  hydrogen 
stoichiometric  ratio,  a  bypass  blower  may  be  required. 

The  relative  humidity  and  the  water  vapor  mass  fraction  at  the 
anode  are  shown  in  Fig.  14.  No  water  vapor  is  condensed  in  the 
anode  flow  channel  due  to  the  high  current.  The  relative  humid¬ 
ity  decreases  with  an  increase  of  the  stack  current,  resulting  in  a 
decrease  of  the  water  vapor  mass  fraction.  Both  of  them  will  affect 
the  ejector  performance,  which  implies  that  fluctuations  in  the 
anode  channel  will  feedback  to  the  fuel  cell  stack  via  the  ejector. 

7.  Conclusions 

A  new  theoretical  model  for  the  convergent  nozzle  ejector  used 
in  the  anode  recirculation  line  of  PEM  fuel  cell  systems  was  pro¬ 
posed  in  this  paper.  A  concave  2D  velocity  function  is  constructed 
by  considering  the  working  conditions  and  flow  characteristic 
of  the  PEM  ejector.  A  method  for  determining  the  exponent  of 
the  velocity  function  nv  was  presented  by  using  the  numerical 
CFD  technique  with  data  regression  and  parameter  identification 
treatments. 

Based  on  the  2D  velocity  function,  the  proposed  modeling 
technique  can  provide  more  accurate  evaluations  of  the  ejector 
performances  than  the  conventional  ID  “constant  pressure  mix¬ 
ing”  and  “constant  area  mixing”  ejector  modeling  theories.  The  new 
model  was  capable  of  analyzing  the  ejector  performances  such  as 
the  primary  mass  flow  rate,  the  secondary  mass  flow  rate,  the  recir¬ 
culation  ratio  and  the  ejector  exit  states  in  not  only  the  critical  mode 
but  also  subcritical  operational  mode. 


6.2.  Fuel  cell  performance 
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Performances  of  the  PEM  fuel  cell  system  are  analyzed  by  vary¬ 
ing  the  stack  currents  as  shown  in  Figs.  12-14.  In  Fig.  12,  the 
performance  properties  of  the  fuel  cell  stack:  voltage,  power  and 
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